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Bad metal properties have motivated a description of the parent iron pnictides as correlated metals 
on the verge of Mott localization. What has been unclear is whether interactions can push these 
and related compounds to the Mott insulating side of the phase diagram. Here we consider the 
iron oxychalcogenides La202Fe20(Se,S)2, which contain an Fe square lattice with an expanded unit 
cell. We show theoretically that they contain enhanced correlation effects through band narrowing 
compared to LaOFeAs, and we provide experimental evidence that they are Mott insulators with 
moderate charge gaps. We also discuss the magnetic properties in terms of a Heisenberg model with 
frustrating J1-J2-J2 exchange interactions on a "doubled" checkerboard lattice. 

PACS numbers: 71.10.Hf, 71.27. +a, 71.55.-i, 75.20.Hr 



Introduction: Iron pnictides are the first non-coppcr- 
based compounds to exhibit high-T c superconductiv- 
ity [J Q, and have therefore attracted considerable in- 
terest. How strong the electron correlations are in the 
parent iron pnictides is a subject of extensive current 
discussion. One viable description places these com- 
pounds at the boundary of itinerancy and Mott local- 
ization. The motivation for this incipient Mott picture 
comes from the observed bad metal properties 0, 3] , and 
from first-principles calculations [jl HI and related 
considerations. Optical conductivity has in the mean- 
time shown a sizable suppression of the Drude weight 
10 12], as well as temperature-induced spectral- weight 



this is the same as in an FeAs layer of LaOFeAs or an 
FeTe layer of the compound FeTe HI, 111- In addition, 



transfer going to energies of the eV range 

m ESQ 

Inelastic neutron scattering experiments have provided 
complementary support. They have not only observed 
zone-boundary high-energy spin waves 15j, but also 
shown 15|, |l6[ that the total (ordered plus fluctuating) 
spin spectral weight is, for instance, about 1.2 /ie/Fe in 
CaFe2As2. Such a large spin spectral weight implies that 
the low-energy spin degrees of freedom come from elec- 
tronic excitations not only close to the Fermi energy but 
also far away from it. All these are defining properties 
of metals on the verge of Mott localization. To establish 
the incipient Mott picture, however, it is important to 
identify the Mott insulating part of the electronic phase 
diagram in the iron pnictides and related compounds. 



Towards this end, we consider the iron oxychalco- 
genides La202Fe20Sc2 and La202Fe20S2. These sys- 
tems are built from stacking layered units La202 and 
Fe20(Se,S)2- We will describe the Se case below, with 
the understanding that the S case is similar unless other- 
wise stated. Fig. Q] (left panel) shows the crystal structure 
with space group I4/mmm as given in Ref. [13] ■ Each 
Fe20Se2 layer contains a square lattice of Fe ions, and 



Fe 2+ is the nominal valence as in the latter materials. 
The Fe20Se2 layer is unique in that the Fe ions in each 
plaquette are alternatively linked by an in-plane oxygen 
ion and by two selenium ions buckled on the two sides 
of the Fe square lattice. The inter-atomic distances are 
rf Fo _ Fc = 2.884 A, dpo-o = 2.039 A, and d Fc -Sc = 2.722 A. 

An important observation is that Fe-square-lattice unit 
cell of La202Fe 2 OSe2 is slightly larger (by about 1%) 
than that of LaOFeAs and considerably larger than that 
of either FeTe (by over 6%) or FeSe (by over 8%) OGlJ. 
This raises the possibility for a narrower bandwidth (oc i) 
and, correspondingly, a larger strength of the normal- 
ized electron correlation, U jt. We have further been mo- 
tivated by the indication of antifcrromagnctic ordering 
(AFM) in the early work of Mayer and co-workers [13] ■ 

Band narrowing: We have determined the elec- 
tronic structure using the WIEN2k linearized augmented 
plane wave method [20| based on density functional 
theory (DFT). A generalized gradient approximation 
(GGA) [21| was used to treat exchange and correlation. 

In Fig. [2 we plot the projected density of states (DOS). 
In both compounds, wc find that the 3d electrons on Fe 
contribute most to the DOS near the Fermi energy. The 
Fe d-electron DOS is mostly confined between —2 eV and 
1.2 eV. This represents considerable narrowing of the Fe 
d-electron bands compared to LaFeAsO, where it occurs 

(There are also some 



22] 



between —2.2 eV and 2 eV 
differences between the two systems in the DFT-derived 
DOS near the Fermi energy, but the effective U/t primar- 
ily depends on the overall 3d-electron bandwidth.) The 
results for La202Fe 2 OS2 [dashed line of Fig. Efa)] are 
similar. A comparable degree of 3d-electron band nar- 
rowing exists when the iron oxychalcogenides are com- 
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FIG. 1: (Color online) (Left panel) Crystal structure 
of La202Fe20Se2; (Right panel) XRD patterns for the 
La203Fe2Se2 and La20sFe2S2 samples. The cell parame- 
ter values are a = b = 4.085 A and c = 18.605 A for 
La 2 3 Fe 2 Se2 and a = b = 4.042 A and c = 17.929 A for 
La20sFe2S2. Atomic positions in the unit cell of a paramag- 
netic state are as follows: La (1/2,1/2,0.18445), Fe (0.5,0,0), 
Se (0,0,0.09669), 01 (1/2,0,1/4), 02 (1/2,1/2,0). 



pared to the iron chalcogenides FeTe and FeSe [231 ] . 

Mott insulating behavior: The narrower Fe 3d elec- 
tronic bands point to enhanced correlation effects. To ex- 
plore this, we have synthesized these materials following a 
similar procedure to that taken in Ref. [171 ] and measured 



their transport and magnetic properties. Polycrystallinc 
samples with nominal composition La203Fe2M2 (M = S 
or Se) were prepared by conventional solid state reaction 
using high purity La2C>3, Fe and S (or Se) powder as 
starting materials. The samples were characterized by 
powder X-ray diffraction (XRD) (right panel of Fig. [T]) 
with Cu Kcv radiation (A = 1.5418 A) at room tempera- 
ture. 

The electrical resistivity as a function of temperature 
is shown in Fig. Ifa). For both M = Se and M = S, the 
behavior is insulating with activation energy gaps [cf. 
inset of Fig. E£a)] of about 0.19 eV and 0.24 eV, respec- 
tively. Figure [HJb) shows the temperature dependence of 
the bulk magnetic susceptibility x{T), which exhibits a 
broad maximum around 120 K. The rapid decrease, at 
93 K and 105 K, for M = Se and M = S respectively, 
is naturally ascribed to antiferromagnetic ordering in the 
Fe d-electron moments. At the lowest measured temper- 
atures, x(T) shows a small increase for La2C>3Fe2Se2 but 
a small decrease for La2C>3Fe2S2. This may be due to dif- 
ferent impurities, whose amounts must be small as there 
is no trace of them in the XRD patterns. 

The results provide evidence that both La203Fe2Se2 
and La2C>3Fe2S2 are Mott insulators, i.e. correlation- 
induced insulators with low energy spin excitations that 
order antifcrromagnctically at low temperatures. 

Magnetism: To investigate the Mott insulating state, 
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FIG. 2: (Color online) The partial DOS for paramagnetic 
La202Fe20Se2 (a) and LaOFeAs (b). Also shown in Panel 
(a) is the Fe 3d partial DOS for La 2 2 Fe20S2. In our DFT 
calculation, the energy threshold separating the localized and 
non-localized electronic states is chosen to be —6 Ryd. The 
muffin-tin radii are: 2.37 ao (Bohr radius) for Lanthanum, 
2.03 ao for iron, 2.40 ao for selenium, and 1.80 ao for oxygen. 
The criterion for the number of plane waves is chosen to be 
i?MT * ^ max = 7 and the number of fc-points is 32 x 32 x 8 
for these paramagnetic-state calculations. 



we extend our DFT calculations to the cases of magneti- 
cally ordered states using the GGA+U method. We con- 
sider seven magnetic ground states, the building blocks 
for which are illustrated in Fig. SJ Each crystalline unit 
cell consists of two layers of Fe20Se2, whose coordinat- 
ing geometry is rotated with respect to each other by 
90° along the c-axis. Correspondingly, the seven spin 
states are formed with the spin alignment as follows: FM 
((a)+(a)), AFM1 ((b)+(b)), AFM2 ((c)+(a)), AFM3 
((d)+(a)), AFM4 ((e)+(e)), AFM5 ((f)+(f)), and AFM6 
((e)+(e)'), where the notation ((spinl)+(spin2)) repre- 
sents only the spin configuration rather than the underly- 
ing atomic coordination within each layer and the symbol 
(c)' represents the spin configuration obtained by rotat- 
ing the configuration (e) in Fig. |U by 90°. The energy 
differences between these spin ordered states for various 
on-site Hubbard interaction U are summarized in Tablefl] 
(None of the magnetically-ordered states at U = has a 
gap in the DOS at the Fermi energy. Within GGA+U, 
an insulating gap does develop for sufficiently large value 
of U, for some of the magnetic ordering patterns.) Based 
on considerations similar to the FeAs case 0, l26l - 28 1. we 
model the spin Hamiltonian in terms of nearest-neighbor 
(n.n.) exchange interaction (Ji) and the next-nearest- 
neighbor (n.n.n.) exchange interaction across an Fe- 
square plaquette with two buckled Se ions (J2) and the 
n.n.n. interaction across an Fc-square plaquette contain- 
ing the O-ion (J^)- Even though the spin system is ex- 
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TABLE I: Relative energies AE (meV/unit cell) of the FM 
(ferromagnetic) and six AFM states of La202Fe20Se2 ob- 
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tained from GGA+17 calculations. 






U (eV) 


Ji/S 2 


J2/S 2 


J2/S 2 


a 









-3.53 


2.10 


9.23 


139.5 






1.5 


2.53 


-2.40 


6.17 


5.5 




y — m = s ^ S \ N 


3.0 


3.77 


-1.11 


4.79 


1.27 


i i 


M = Se 




4.5 


3.38 


-0.78 


3.28 


0.63 



50 100 150 200 250 300 

7(K) 



FIG. 3: (Color online) (a) The electrical resistivity (p) vs. 
temperature (T) for La20aFe2S2 and La20sFe2Se2. Inset: 
lnp vs. 1/T; (b) The temperature dependence of susceptibil- 
ity, \, measured at 3 Tesla. The three open points, measured 
at 2 K, 50 K, and 150 K, respectively, are determined from 
the slope of the M(H) curves. 
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FIG. 4: (Color online) Schematic representation of the 
spin building blocks for the magnetically ordered states of 
La202Fe20Se2. To be consistent with the effective model, 
only Fe atoms are shown. The arrows represent the spin ori- 
entation. The black solid, green dashed, and blue solid lines 
stand for Ji, J2, and J2 spin exchange paths. 



pected to be largely Heisenbcrg-type, the exchange pa- 
rameters can be extracted by fitting the ground-state en- 
ergies in terms of the Ising counterpart [HI [HI . We sum- 
marize the extracted exchange interactions in Table [TT1 

We now turn to the effective J \- J 2- J 2 Hcisenberg 
model defined on the two-dimensional doubled checker- 
board lattice shown in Fig. [5ja). The classical Hamilto- 
nian is H = Jij^i ■ Sj, where Si is a 3D vector 



TABLE II: Values of spin exchange parameters J1/S 2 , J2/S 2 , 
J2/S 2 (meV) for La20 2 Fe20Se2 estimated from GGA+U cal- 
culations. Positive (negative) J's correspond to antiferromag- 
netic (ferromagnetic) exchange. The last column lists the 
"error bar" in energy (meV/unit cell): a = [J2iLi( AE ~ 



6 and A£« and AE<$ 



AE^f/iN - l)] 1 ' 2 , where N 
the relative ground-state energies in DFT and in the extracted 
effective spin models, respectively, a is large for (7 = 0, but 
becomes relatively small for moderate values of U . 



with magnitude S. To study the T = phase diagram, 
for J 1 > 0, we rewrite the Hamiltonian via a summa- 
tion over two sets of inter-penetrating plaquettes shown 
in Fig. E^a): 



H = -{J2 + J' 2 )NS' 2 
+ J 2m [(S! + S 3 ) 2 



(S 



m=A,B Orn 
2 + Sa^ 2 



(1) 



where J\a 



aJi, Jib = (1 - a)Ji; J 2 a(b) = ^2(^2) 



Ji a(b)- Here ^ a 1 is a dummy control param- 
eter [29(, and N refers to the total number of plaque- 
ttes. The phase diagram is shown in Fig.[5jb). There are 
the Neel ordered AFM [for < a < 1, J 2 < aJi, J' 2 < 
(1— a) Ji] and the collinear (ir, 0)/(0, n) ordered AFM [for 
< a < 1, J2 > aJi, J' 2 > (1 — a)Ji] states. Unlike the 
square lattice J1-J2 model, there also exists a plaquette 
(7r/2,7r/2) AFM phase, which is stable for ^(J^) > Ji 
and J' 2 {J2) < 0. The transition between any two of the 
four ordered regions in the phase diagram is first order. 
Interestingly, we find that along the red solid boundaries 
in Fig. [5{b) , the model loses any long-range magnetic or- 
der since the ground-state is infinitely degenerate. Along 
this boundary, J2 = J' 2 = 1/2 corresponds to the maxi- 
mally frustrated point of the J1-J2 model, and ^2(^2) = 
corresponds to the checkerboard AFM model. The phase 
diagram for J\ < can also be obtained by rotating all 
the spins in one sublattice of the underlying square lat- 
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tice (see Fig. E^a)) S, — > — S% and sending Ji — > — J\. It 
is similar to the one for J\ > 0. The only difference is 
that the Neel ordered (it, n) phase is replaced by the FM 
phase ordered at (0,0). 

The magnetic ground states expected from the ex- 
change couplings given in Table [TT| are marked in 
Fig. EKb). For U = 1.5 eV and U = 3.0 eV, the ground 
state is the (tt/2,tt/2) plaquette AFM state, while for 
larger U, i.e., U = 4.5 eV, the standard (n,ir) Neel state 
is energetically favorable. Taking into account the stack- 
ing structure of the material, we see that these results arc 
consistent with the DFT energies listed in Table Q] the 
AFM6 state has the lowest energy for U = 1.5 eV and 
U = 3.0 eV, and the AFM1 state has the lowest energy 
for U = 4.5 eV. 

To summarize, we studied the iron oxychalcogenidcs 
La202Fc20(Se,S)2, whose Fe-square lattice is expanded 
compared to that of the usually considered iron pnictides 
and chalcogcnides. We have theoretically demonstrated 
that their Fe 3<i-electron bands are narrower than the 
usual cases. The corresponding enhancement of correla- 
tion effects promotes the Mott insulating state, and this 
is demonstrated experimentally. Our demonstration of a 
Mott insulating phase in systems closely related to the 
iron pnictides extends the electronic behavior of these 
systems to a new regime. In addition, our results support 
the notion that the iron pnictides/chalcogenides possess 
intermediately strong electron correlations and are not 
too far away from Mott localization. 
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FIG. 5: (Color online) (a): The doubled checkerboard lattice 
for the Hamiltonian defined in Eq. (TjQ). It consists of two sets 
of inter-penetrating plaquettes, indicated as blue (A) and red 
(B) dashed squares. It can be decomposed into two (red and 
black circle) sublattices. (b): The T = phase diagram. 
The occupied symbols mark the ground states with exchange 
couplings listed in Table [TTt with U — 1.5 eV (triangle); U = 
3.0 eV (square); U = 4.5 eV (circle). 
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